"Weak Quantum Chaos" and its resistor network modeling 
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Weakly chaotic or weakly interacting systems have a wide regime where the common random 
matrix theory modeling does not apply. As an example we consider cold atoms in a nearly integrable 
optical billiard with displaceable wall ("piston"). The motion is completely chaotic but with small 
Lyapunov exponent. The Hamiltonian matrix does not look like one taken from a Gaussian ensemble, 
but rather it is very sparse and textured. This can be characterized by parameters s and g that 
reflect the percentage of large elements, and their connectivity, respectively. For g we use a resistor 
network calculation that has a direct relation to the semi-linear response characteristics of the 
system, hence leading to a novel prediction regarding the EAR of cold atoms in optical billiards 
with vibrating walls. 



I. INTRODUCTION 

So called "quantum chaos" is the study of quantized 
chaotic systems. Assuming that the classical dynamics 
is fully chaotic, as in the case of a billiard with convex 
walls (Fig.[T]), one expects the Hamiltonian to be like a 
random matrix with elements that have a Gaussian dis- 
tribution. This is, of course, a sloppy statement, since 
any Hamiltonian is diagonal in some basis. The more 
precise statement behind random matrix theory (RMT) 
is following [iHll ■ Assume that there is a Hamiltonian % 
that generates chaotic dynamics, and consider an observ- 
able F that has some classical correlation function C(i), 
with some correlation time t^ . Then the matrix represen- 
tation Fnm in the basis of % looks like a random banded 
matrix. The bandwidth is ?ipianck/^R- If is small, such 
that the bandwidth is large compared with the energy 
window of interest, then the matrix looks like taken from 
a Gaussian ensemble. 

Our objective is to analyze the energy absorption rate 
(EAR) of billiards with vibrating walls, which is related 
to past studies of Nuclear friction [3-|9j. However, our 
interest is focused on 2D optical billiards [lol - [l^ whose 
geometrical shape can be engineered. In this problem Ti. 
is the Hamiltonian of the non-driven billiard, while Fnm 
is the perturbation matrix due to the wall displacement. 
If the driving is not too strong we expect a linear relation 
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where / is the RMS value of the vibrating wall velocity. If 
one further assumes that the billiard is strongly chaotic, 
then G = Go can be determined (Eg. (|M)) ) from simple 
kinetic considerations as in 04^, leading to a variation 
of the so-called Wall-formula. Note that there is a strict 
analogy here with the Drude formula and the Joule law. 

We consider completely chaotic billiards [3, , with 
no mixed phase space, but we assume that they are only 
weakly chaotic [17h24| . This means that is much larger 
than the ballistic time ti^. Consequently, the EAR coef- 
ficient is 



with g 1. In the classical analysis g = gc is related 
to classical correlations between the collisions with the 
vibrating walls. In the quantum analysis the first ten- 
dency is to assume g ~ gc- In contrast to that wc would 
like to highlight the possibility to observe g ^ gc- This 
is the case if we have weak quantum chaos (WQC) cir- 
cumstances, in which the traditional RMT modeling does 
not apply, meaning that i^„m does not look like a typical 
random matrix. Rather, the distribution of its elements 
is log- wide (resembles a log- normal distribution), and it 
looks very sparse, as expected from Refs.[25l-l28j. Conse- 
quently, the analysis of the EAR has to go beyond the 
familiar framework of linear response theory (LRT). 

WQC circumstances are encountered in the analysis of 
any weakly chaotic or weakly interacting system. In the 
WQC regime the matrix Fnm is formed of elements that 
have a log-wide distribution. The implied sparsity is im- 
portant for the analysis of the EAR [3: l23l , a s expected 



from semi-linear response theory (SLRT) |29l - l3l| . The 
main idea behind the theory is demonstrated in Fig.[2] 
one observes that the energy absorption process requires 
connected sequences of transitions between the energy 
levels of the system. Accordingly, the calculation of EAR 
requires a semi-linear resistor network calculation. 

Wc can characterize the sparsity of the perturbation 
matrix by parameters s and gs that rcfiect the percent- 
age of large elements, and their connectivity, respectively. 
The parameter gs is defined through a resistor network 
calculation, and has a direct relation to the semi-linear 
response characteristics of the system, namely 
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For a strictly uniform matrix g^ = s = 1, for a Gaus- 
sian matrix s = 1/3 and ^ 1, while for sparse matrix 
s,gs <C 1. We would like to explore the dependence of 
gc and gs on the parameters u and h of the system. 
Disregarding the physical motivation, this exploration is 
mathematically interesting, because it introduces a "re- 
sistor network" perspective into RMT studies. Hence, it 
is complementary to the traditional spectral and inten- 
sity statistics investigations. 
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A. Generic parameters 

For a billiard of linear size L that has walls with radius 
of curvature i?, as in FiglU the Lyapunov time and the 
ballistic time are 

tn = R/vj^ [Lyapunov time] (4) 
tL = L/v-E, [Ballistic time] (5) 

where = {2E/mY^'^ is the velocity of the particle. The 
quantization introduces an additional length scale into 
the problem, the de Broglie wavelength 
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Accordingly, the minimal model for the purpose of our 
study is featured by two small dimensionless parameters: 

u= L/R = [degree of deformation] (7) 
h = Ae/L = [scaled Planck] = 2tt /{k^L) (8) 

With the two classical time scales and one may 
associate two frequencies, while quantum mechanics adds 
an additional frequency that corresponds to the mean 
level spacing: 



Ao = 27r/tH = {h/2TTY-^ Al 
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where d=2 is the dimensionality of the billiard. The 
WQC circumstances that we would like to consider are 
characterized by the following separation of scales: 



Ao, Ar < Al 
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However, this is not a sufficient condition to observe 
WQC. The identification of the WQC regime in the (u, h) 
space is an issue that we would have to address. 



B. Detailed outline 



The model.— In Sec.(|TT| we define the model and ex- 
plain the numerical procedure. Schematically the Hamil- 
tonian of the system can be written as 
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where "Ho describes the undeformed rectangular box, and 
U describes the deformation of the fixed walls, and F is 
the perturbation due to the displacement fit) of the mov- 
ing wall (piston) . The geometry of the billiard is charac- 
terized by u, while h defines via Eq.® the energy window 
of interest. 



Eigenstates.— Given u and h we find the ordered 
eigenenergics £"„ of the Hamiltonian "H, within the en- 
ergy window of interest. This is done using the bound- 
ary element method [32|. A representative eigenstate is 
presented in Fig[3l If the deformation is small it is mean- 
ingful to represent it in the basis which is defined by Hq. 
See for example Fig21 As the deformation becomes larger 
more and more levels are mixed as demonstrated in Fig[S] 

Perturbation matrix.— Once wc have the eigen- 
states we calculate the matrix elements Fnm of the per- 
turbation term _F, within the energy window of interest. 
An image of a representative matrix is shown in Fig.|6l 
The distribution of its elements is definitely not Gaus- 
sian, as shown in Fig[7l In fact we see that the statistics 
of \Fnm\'^ resembles a log- normal distribution. 

Bandprofile.— In order to characterize the bandpro- 
filc of the perturbation matrix we define in Sec. pi)) spec- 
tral functions Ca [r) and {r) that are displayed in Fig[8l 
We explain how Ca is semiclassically related to the power 
spectrum C{u)) of the collisions, and how Cs gives an in- 
dication for the sparsity of the matrix. 

Sparsity.— Wc characterize the sparsity and the tex- 
ture of the matrix by a parameter gs that is defined in 
Sec. (|IIip . The numerical results for are presented in 
Fig-El This characterization cares about the connectiv- 
ity of the matrix elements, and is based on a calcula- 
tion of a resistor network average. Some further details 
with respect to the resistor network average are given in 
App.®. 

Classical analysis.— In Sec. (jIV[) we provide a de- 
tailed analysis of the classical power spectrum C{uj). 
In particular, we derive an expression for the zero fre- 
quency peak, which reflects the long time correlations of 
the bouncing trajectories in our weakly chaotic billiard. 
App. (|B|) provides optional perspectives with regard to 
this classical calculation. 

Quantum analysis.— In Sec. dV]) we use 1st order 
perturbation theory for the analysis of the eigenstates of 
the deformed billiard, and hence get an approximation for 
C(w). The validity of this approximation is very limited. 
We therefore fuse perturbation theory with semiclassical 
considerations in Sec. (|VI|) . This allows to obtain some 
practical approximations for g^ and gs- 

The WQC regime.— Eventually we turn to define 
the borders of the WQC regime in the (u, h) parameter 
space. This is also an opportunity to make a connection 
with previous works that concern spectral and intensity 
statistics for such type of weakly chaotic billiards. 

Implications.— The relevance of the resistor network 
analysis to the EAR calculation is clarified in Sec. ljVIIip . 
and the experimental feasibility of observing the implied 
SLRT anomaly is discussed in Sec. ljIXp . A broader per- 
spective with respect to EAR predictions is presented 
in Sec.(jX|. In particular, we clarify how to bridge be- 
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twcen what looks like contradicting results with regard 
to diffusive and difFractive systems, as opposed to ballis- 
tic billiards. 



II. THE MODEL, NUMERICS 



The Hamiltonian of the system is 

^to..i = 1^ + W,^^{x,y) + Wf{x-f{t)Ds{y)) (15) 

We write it formally as in Eg. (fT3l) . where T-Lq describes 
the undeformed rectangular box x Ly. The ballistic 
time is defined as tj^ = Lx/v^, which is the typical time 
for successive collisions with the piston. The term U in 
the Hamiltonian is due to a deformation Du{y) of the 
left static wall. The amplitude of this deformation is 
Do ~ Ly/R, while u = Ly/R is conveniently defined as 
the dimensionless deformation parameter. The driving 
term ~f{t)F is due to the deformation f{t)Df{y) of the 
right wall, leading to the identification 



The associated power spectrum is: 



F = Dfiy)W'fix) 



(16) 



For parallel displacement of a "piston" we set Df{y)=l. 
Note that f{t)Df unlike £>„ is assumed to be small com- 
pared with the de Broglie wavelength. 

For a billiard system with "hard" walls the potential 
is zero inside the box, and become very large outside of 
the box. Accordingly, it is assumed that Wf{r) = for 
r < 0, with a steep rise as r becomes > 0. Accordingly, 
the penetration distance upon collision is much smaller 
compared with the linear dimension L of the box. Say 
that the force which is exerted on the particle by the 
piston is Wj:{r) = Fq for r > and zero otherwise, then 
it is assumed that E/Fq <^ L, where E is the kinetic 
energy of the particle inside the box. Below we take the 
limit _Fo — oo. 

Following [ssl - lssj we discuss the definition and the cal- 
culation of the spectral function C{uj) that describes the 
fiuctuations of F in the non-driven billiard system. We 
first discuss the classical and then turn to the quantum. 

Classical.— In the classical context the Hamilto- 
nian y. can be used to generate a trajectory {tj ,yj,9j), 
where j labels successive collisions with the piston at 
{0,yj) with incident angle 9j. Consequently, the asso- 
ciated F{t) consists of impulses of height Fq whose du- 
ration is 2mvE cos{dj)/ Fq. In the hard wall limit one can 
write formally 

F{t) = ^2m^;, cos((?,) %(2/,) J(t-tj) (17) 
j 

Assuming ergodic motion, the auto-correlation function 
of F{t) can be calculated from the time dependence of a 
single trajectory that has some long duration 
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If we regard the impulses as uncorrelated we get the result 
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which holds in the uj oo limit. More details about 
this calculation and its refinement will be presented in 
SecdrvT). 



Quantum.— The unperturbed energy levels of the 
rectangular box are 
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with the mean level spacing 
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For a given deformation we diagonalize and find the 
ordered eigenenergies £"„ with n = 1,2,3,..., within an 
energy window of interest which is characterized by the 
dimensionless parameter h. This is done using the bound- 
ary element method [32|. Each eigenstates '0(a;, y) is rep- 
resented by a boundary function ^p[y) = ip'(0,y), where 
the normal derivative is with respect to x at the position 
a;=0 of the piston. Consequently, the matrix elements 
of F are 
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— / v'^")(y)^(")(2/) Df{y)dy (24) 



Given i^„,„ one can calculate the quantum mechanical 
version of the spectral function 
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where it is implicit that the delta functions have a finite 
smearing width related to the measurement time too ■ and 
an average over the reference state {En ^ E) is required 
to reflect the associated uncertainty in energy. 

Correspondence.- For a chaotic system, if the cor- 
relation time is short, one expects quantum-to-classical 
correspondence (QCC) with regard to C{u!). It follows 
from Eq. ([25]) that this spectral function should reflect 
the bandprofile of the perturbation matrix Let us 

express this observation in a convenient way that allows 
a practical procedure for numerical verification. We cal- 
culate Fnm in an energy window of interest, and define 
the associated matrix 



C{t) = {F{t)F{0)) = F{t)FiO) 



(18) 
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The bandprofile Ca {r) is defined by the average of the el- 
ements Xnm along the diagonals n—m = r. In the same 
way wc also define a median based bandprofile Cs{r). 
Given that the mean level spaeing Aq is small compared 
with the energy range of interest, the correspondence be- 
tween Cq,„(a;) and Cci(w) can be expressed as: 



Ca{n - m) 
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In particular, it follows from Eq.([2T|) that the unrestricted 
average value of the elements Xnm is 
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In fact this result can be established without relaying on 
QCC considerations via a sum rule that we discuss in 
Sec. (|VI|) . and the same result is also obtained from the 
zero order evaluation of matrix as described in Sec. (|V| . 
Whenever applicable we re-scale the numerical results 
with respect to this reference value. 

The applicability of the QCC relation Eq. ((27| . to the 
analysis of our billiard system is confirmed in Fig[8l down 
to very small frequencies. We also see that 



Csir) < Cair) 



(29) 



where the value on the right is obtained in the limit 
r — > oo . The inequality Cs ^ Ca means that the value of 
the typical matrix element is very small compared with 
the average value. We are therefore motivated to define 
notions of sparsity and texture in Sec. pilj) . 



III. SPARSITY AND TEXTURE 

For strongly chaotic systems the elements within the 
band have approximately a Gaussian distribution. But 
for WQC the matrix becomes sparse and textured as 
demonstrated in Fig.Hl These features go beyond the 
semiclassical analysis of the bandprofile. The sparsity is 
related to the size distribution of the in-band elements: 
Loosely speaking one may say that only a small fraction 
(s) of elements are large, while most of the elements are 
very small (for a precise definition of s see below). The 
texture refers to the non-random arrangement of the mi- 
nority of large elements. 

In the WQC regime the size distribution of the in- 
band elements becomes log-wide (approximately log- 
normal) as seen in Fig.[71 This is reflected by having 
Cs{r) ^ Ca{r) as seen in Fig[8] Accordingly, an optional 
measure for sparsity is the parameter q which is defined 
as the ratio of the median to the average. 

The sparsity and the texture of are important for 
the analysis of the energy absorption rate as implied by 
SLRT. Accordingly, it is physically motivated to char- 
acterize the sparsity by a resistor network measure gs 
that reflects the connectivity of the elements, and hence 



has a direct relation to the semi-linear response char- 
acteristics of the system. The precise deflnitions of gs is 
given below. For a strictly uniform matrix gs = s ~ I, for 
a Gaussian matrix s = 1/3 and (7s ~ 1, while for sparse 
matrix s,gs <SC 1. The dependence of the sparsity on the 
energy and on the degree of deformation is demonstrated 
in Fig[Sl and is related to the mixing of the levels in Fig[S] 

Definition of s.— Deflne a matrix X whose elements 
are Xnm — l^nmP- Associate with it an untextured ma- 
trix X"'" and a uniformized matrix X""' that have the 
same bandproflle C'(r). The former is obtained by per- 
forming random permutations of the elements along the 
diagonals, while the latter is obtained by replacing each 
of the elements of a given diagonal by their average. The 
participation number (PN) of a set {Xi} is deflned as 
(^^ A"i ) Xf ; and reflects the number of the large el- 

ements. Here the index is i = {n,m). The PN of {Xnm} 
counts the number of large elements in the matrix. The 
PN of {Xjj^n} counts the number of the in-band elements. 
Accordingly, the ratio constitutes a measure for sparsity: 



s[X] 



PN [X] 

PN 



(30) 



It should be clear that X and X""" have the same s 
but only the former might have texture. So the next 
question is how to deflne texture avoiding a subjective 
visual inspection. 

Definition of g^.— Complementary information 
about the sparsity of the matrix, that takes into account 
the texture as well, is provided by the resistor network 
measure gs- Coming back to A"„ 
it a matrix 



we can associate with 



Gnm = 2F(n-m) 



Xn 



(n — m)^ 



(31) 



where -P"(^) = 1 is a weight function whose width 
should be quantum mechanically large (i.e. ^ 1) but 
semiclassically small (i.e. < the bandwidth). With such 
choice the Gnm are proportional to the Fermi-golden rule 
transition rates that would be induced by a low-frequency 
driving —f{t)F. Optionally we can regard these Gnm as 
representing connectors in a resistor network, as in Fig[2] 
The inverse resistivity of the strip can be calculated using 
standard procedure, as in electrical engineering, and the 
result we call {{X))s. For more details see App.(jX|. It is 
useful to notice that if all the elements of X are identi- 
cal then {{X))s equals the same number. More generally 
{{X))s is smaller than the conventional algebraic aver- 
age {{X))a (calculated with the same weight function). 
Accordingly, the resistor network quantity {{X))s can be 
regarded as a smart average over the elements of X , that 
takes their connectivity into account. Consequently, it is 
natural to define a physically motivated resistor-network 
measure for sparsity and texture: 
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(32) 
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One can show that {{X))s is strictly bounded from be- 
low by the harmonic average. In practice the geometric 
average or the median provide better lower bounds. In 
the RMT context a realistic estimate for {{X))s can be 
obtained using a generalized variable-range-hopping pro- 
cedure (see [3l| for details). 



Additional definitions.— If the elements Xnm have 
a well defined average {{X))oo in the limit of infinite 
truncation, then it is convenient to define 



_ {{X))a 
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(33) 
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Later, in Sec. ljVIIip . we discuss the physical significance 
of (7, and identify it as the dimensionless absorption coef- 
ficient. In particular, is identified as the dimensionless 
absorption coefficient in the classical calculation, which is 
determined by taking into account classical correlations. 
In the quantum case wc have an additional suppression 
factor gs due to the sparsity of the perturbation matrix. 



IV. CLASSICAL ANALYSIS OF C{uj) 

Let us assume that we have a collision in angle 9 with 
the flat piston. The force which is exerted on the particle 
during the collision is Fq, such that the impact is 



2mwE cos(6') 



(35) 



Consequently, the force F{t) looks like a train of spikes 
as in Eq. (fT7|) . Note that the duration of a collision is 
qe/Fo. In the absence of deformation the time distance 
between the spikes is 



2L, 



Te = 



(36) 



where 9 is constant of the motion. For the following 
calculations it is useful to define the following averages: 
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(39) 
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If we have a very small u the effect would be to er- 
godize 9 with some rate 'yg. After time t the number 
(#) of collisions is t/rg and consequently the deviation of 
the perturbed trajectory is multiplied by (1 -I- L/R)"^ m 



exp ((uE^/i?) cos(0)). Accordingly, the instability expo- 
nent is 



70 



70 



R 
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(41) 



For sake of generality we have added a background term 
7o. This background term would arise if the upper or 
lower walls were deformed, or if the potential floor were 
not flat. A non zero 70 is unavoidable in a realistic sys- 
tem. As we shall see shortly the effect of the deformation 
is twofold. The primary effect is to ergodize 9, and the 
secondary effect is to modify the small uj spectral content 
of the fluctuations. 

Power spectrum.— Let us deflnc F0{t) as the tem- 
poral "signal" which is associated with a trajectory that 
starts at the piston with 9 collision angle. This signal 
consists of delta spikes, the first one being qg5{t). The 
correlation function can be expressed as 



{F{0)Fit))^{:i^Fo{t) 



correlations 



(42) 
(43) 



where the first term represents the self-correlation of the 
spikes. It is convenient to subtract from {F{0)F{t)) its 
global offset, and to define the correlation function as 

Cit) = {FiO)Fit)) ~ {Fr (44) 

The associated power spectrum is the Fourier transform: 



(45) 



^ 2n5{uj) 

TO 



Note that Fq(uj) is the FT of Fg{t). It is not the same as 
F{ll>) of Eq. ([T9| . The latter has a random phase due to 
a time displacement of the time origin, while in the case 
of F0{t) the time origin is fixed by the presence of F{0) 
in Eq. (jt2l) . 



Infinite frequency limit.— The first obvious obser- 
vation is that that for large frequencies the power spec- 
trum becomes flat and reaches a constant value that re- 
flects the self-correlation peak of Ea. (|^5|) 



C{oo) = 



(46) 



where Coo is given by Eq. ([57]) . This result, if it is a 
applied to flnite frequencies, is termed in the literature 
"the white noise approximation" . 

Zero deformation.— Let us consider the non- 
deformed integrable billiard. Then the bouncing tra- 
jectories consist of equal spikes and may have an ar- 
bitrary long periods Tg. The Fourier transform of 
Fg{t) — qgS{t—tj), is a reciprocal comb, namely 



Fg{uj) = —S (uj - —n 



(47) 
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The power spectrum is obtained using Eq. pSj) . It consists 
of two components. One component is the zero frequency 
peak which reflects the dispersion of the impact pulses 



Ciuj ~ 0) 



Var 



27r(5(w) 



(48) 



This zero frequency peak would be broadened if the de- 
formation were non zero, as discussed in the next para- 
graph. The second component of the power spectrum 
consists of ballistic peaks at a;„ ~ [tiv^/ Lx)n, that merge 
to Coo in the infinite frequency limit: 



C{uj > 0) 



2n ^1 _ (,,/c,„)2 



(49) 



Fig[5] presents the numerical data for a slightly deformed 
billiard. Disregarding the broadened zero-frequency 
peak, the above zero-deformation result provides a prac- 
tical overall approximation. 

Small deformation.— For small deformation the 
main effect is the broadening of the delta function in 
Eq. (j48|) . Assuming a 6 independent 7. the 5{ijj) is re- 
placed by the Lorentzian (l/7r)7/(cj^ +7^)- Hence, we 
get for small frequencies 



C(w « AO 
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(50) 



We further illuminate this result using a time dependent 
and number variance approaches in App.([B]). If 7 = 70 
is well defined there is a well defined limiting value as 
w — >■ 0. With the identification 7 ~ l/tn one should re- 
alize that the power spectrum at zero frequency is en- 
hanced by factor t^/tj^, hence 



C{uj = 0) 



1 



(51) 



But if 7 is given by Eq. (j41[) we have to perform an ergodic 
average over 6. This becomes interesting if 70 is very 
small or zero, as discussed in the next paragraph. 

Bouncing effect.— It has been proven [l^, [l^ that 
in strictly hyperbolic billiards the time correlation func- 
tion exhibits an exponential decay rate. But if there are 
bouncing trajectories, that do not collide with the de- 
formed surfaces, and can be of arbitrarily long length, 
then a power law decay shows up as in the hard-sphere 
gas and in the Stadium [l^. Our billiard, with 7 as 
given by Eq. (|4T|) . can be regarded as a related variation 
on this theme. Assuming that 70 is very small, the tra- 
jectory has a very long bouncing period when 9 ~ tt/2. 
Consequently, the ergodic average over I/7 generates a 
logarithm factor log(l/7o), or at finite frequency it be- 
comes log(l/a;). Let us be more precise in the 70 = 



case. Averaging over the Lorentzian we get 
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For small lu the above expression can be further simpli- 
fied: 



2Z! 



(53) 



In the zero frequency limit, if 70 is finite, the logarithmic 
factor in the baove expression is replaced by ln[2/(7otR)]. 
Consequently, the result gc = l/u which is implied by 
Ea. ([5T|) is replaced by 



9c = In 



70 



1 

u 



[for 



0] 



(54) 



In the quantum case, that we discuss later, the finite 
level spacing provides an additional lower cutoff Ao that 
"competes" with 70 as discussed in Sec.(|VI|). 



V. PERTURBATION THEORY ANALYSIS 

In this section we shall see what comes out for the 
matrix elements Fnm within the framework of quantum 
perturbation theory to leading order: zero order evalua- 
tion for the "large" elements, and first order perturbation 
theory (FOPT) for the "small" elements. In Sec. (jVI[) we 
shall try to reconcile the perturbation theory results with 
the classical results of Scc. pV[) . 

The small parameter in the perturbative treatment 
is u. The eigenstates n = {n,j:,ny) of the non-deformed 
billiard are 



sin ( Tlx— x\ sin ( ny—y \ (55) 



The deformation profile is 

D^{y) = - (y-e)2 - ^ - [Ly-ef (56) 

In the FOPT treatment the perturbation term in the 
Hamiltonian is calculated using an expression analogous 
to Eq. ([24| , with £)„ replacing D / along the left wall, 
leading to 

(57) 



iL3 



where 



D, = ^ " Duiy) cos (^i^^y ] dy (58) 
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In the numerical analysis we calculate Di^ and hence Unm 
numerically. But here, for presentation purpose, we in- 
troduce a practical approximation: 



T\LIJ 1 



(59) 



In this expression an exponent a=l would arise due to 
the discontinuity of Du{y) at y = 0. However, the effec- 
tive value of a is larger because this discontinuity is very 
small and hardly expressed numerically. Furthermore, 
we would not like to restrict the analysis to the specific 
deformation that had been assumed in the numerics. We 
therefore regard a, for the sake of further discussion, as 
a fitting parameter. 

We can regard the deformation U as inducing scat- 
tering between the Uy modes of the rectangular "waveg- 
uide" . If the box is not deformed (w=0), which is like "no 
scattering", then Uy is a good quantum number. Oth- 
erwise, for non-zero deformation, the levels are mixed. 
The FOPT overlap between perturbed and unperturbed 
states is 



{m\n) 



Um.r}. 



(60) 



Note that by adiabatic continuation we assume in this 
expression an association of perturbed states n with un- 
perturbed states fi={ny,ny). This association holds 
for those levels that are not mixed non-perturbatively. 
Later we discuss the coexistence of perturbative and non- 
perturbative mixing. 

Zero order elements.— We turn to look at Fnm- For 
zero deformation it is block-diagonal with respect to Uy. 
Namely, 



(61) 



mi?. 



Most of the matrix elements are zero, while a small frac- 
tion arc finite. Considering the elements within an energy 
shell setting |7t,| ^ |m| ~ k^^^L, the size of the large el- 
ements is 



IF, 



1 



nrn |0 



li3 



(62) 



In App.([C]) we show that the fraction of elements that 
have this large value is 



TTk^Ly 



(63) 



Consequently, the average value ((|F„,„p))oo of the ele- 
ments is po X |i^nm|o- I^i the more careful calculation of 
App.jC]) we show that 



{{\Fnrn\'))c 



(64) 



in consistency with the semiclassical relation Eq. ([?7)) . 



FOPT elements.— For small u the large size matrix 
elements of F^m are hardly affected by the mixing. But 
at the same time the deformation gives rise to in-band 
small size matrix elements, that would have been zero if u 
were zero. Within FOPT the following approximation 
applies: 



E 

n'm' 



'n\n')Fn'fn.'{'m'\in) 



« Fftffi + (n|m) Fnn + {m\n) Ff, 

— Ffirn (?l|7Tz) {^Ff^fi irrirn) 

Hence, the emerging small elements are 



En 



Effi 



\F^ 



F„ 



Do 



mL* {\n\^ - \m\^) {1 + \ny - my\°') 



(65) 

(66) 
(67) 



(68) 
(69) 



where for simplicity we had assumed ^ Ly ^ L such 
that En ~ 7r2|n|V(2mi2) with \n\ = {nl -f n2)i/2. 

Given an energy window around E, we would like to es- 
timate the typical size of the elements Fnm that connect 
energy levels that have the separation \En — Em\ = OJ- 
Our interest is in small frequencies Aq <C w <C Al. 
Setting ||n| — |m|| ^ w/Al, and Dq « L? / R, and 
\nx—mx\ ^ \ny—my \ ~ \n\ ~ |m| ^ k^L, we get for the 
majority of elements the estimate 



\F 



nm FOPT 



A^ 



1 



mL3 



(k^L) 



3-a 



(70) 



This should be contrasted with the zero order value 
Eq. (|62|) of the large but rare elements: it is much smaller 
whenever the FOPT estimate applies. 



VI. QUANTUM ANALYSIS OF C{uj) 

The QCC relation Eq.(l27]) implies that C^,^{lj) reflects 
the algebraic average over the elements of the matrix 
{l^nmP} along the diagonal \En — E,n\ ^ i^. Our nu- 
merics show that we can trust Eq. ([27]) up to the very 
small frequency Aq. This statement is based on some 
assumptions that should be clarified. 

First we would like to emphasize that both classically 
and quantum mechanically 



C,„(CJ > Al) 



(71) 



In the classical context this value merely reflects the self 
correlation of the spikes of which F{t) consists, and hence 
it is proportional to the ratio between the area (length) 
of the piston and the volume (area) of the box (billiard) . 
In the quantum context it reflects the associated assump- 
tion that well separated eigenstates look like uncorrelated 
random waves, and hence \Fnm\'^ is determined by the 
same ratio as in the classical case. For more details see 
appendices of fssj . 
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QCC condition.— As we go to smaller frequencies, 
correlations on larger time scales become important, and 
the validity of the QCC relation Eq. (P7)) becomes less 
obvious. Recall that due to the bouncing 



a,(c<j-o) 



1 



-Co 



(72) 



Recall also that the matrix elements are strictly bounded 
from above. The maximal value is in fact given by Ea. (|62p 
and accordingly 



n 



(73) 



This has an immediate implication: QCC cannot hold 
globally unless h < u. This requirement can be illumi- 
nated from an optional perspective. The zero frequency 
peak of (7^1 (w) has a width Ap,. This peak cannot be 
resolved by C^^{uj) unless > Ag. Again we get the 
same necessary condition 



h < u, 



[QCC requirement] 



Sum rule.— Extending the discussion with regard 
to Eq. (|7T|) , it is important to realize that the integral 
over Cai{i^) equals Var(F), and accordingly it does not 
depend on u, but only on the ratio between the area 
(length) of the piston and the volume (area) of the box 
(billiard). Note that the height of the zero frequency peak 
is proportional to 1 /u, while its width is proportional to u 
in consistency with this observation. 

In complete analogy, in the quantum analysis the sum 
J2m l-^nmP docs uot depend on u. If the diagonal ele- 
ments can be neglected it follows that (^^^(cj) does not 
depend on u. But if h > u, the zero frequency peak 
cannot be resolved, and the deficiency can be attributed 
to the diagonal elements, in consistency with the FOPT 
analysis. 

FOPT.— In the regime u < ?i it is instructive to 
contrast the lower bound FOPT result which is implied 
by Eq.^U^, with the SC result which is implied by Eq. ((50)) 



C. 



qm-FOPT 



C'qn,-Sc(w) 



3-2q 



1 

h 

A^Ah 
C.2 + A2 



(75) 
(76) 



The lower cutoff Aq in the FOPT expression has been en- 
tered by hand to indicate its existence. It is implicit here 
that the frequency range of interest is A^ ^ w <C Al. 
In the worst case of having a deformation with discon- 
tinuity (a = 1), the ratio between these two results, in 
the frequency range of interest, is as one could expect 
{u/h) <Si 1. We shall discuss the relevance of the FOPT 
and semiclassical expressions below, and also in Sec. (|VIip . 

Evaluation of g^.- Coming back to the regime h < u, 
assuming that the QCC relation Ea. (P7|) can be trusted. 



we deduce that the unrestricted average value of the ma- 
trix elements at energy E is 



{{\Fnn.\'))c 



2tt 



Co 



(77) 



Our interest is in the response characteristics of the sys- 
tem for low frequency driving, which we further discuss 
later in Sec. (|VIIip . We assume that the spectral con- 
tent of the driving is characterized by a cutoff frequency 
Afj < utc < Al . Therefore we look on the band-averaged 
value: 

((|F„„,n), ^ (^) ' - rCU^)dco (78) 
= 9c^ ((|F„„|2))^ (79) 

If QCC holds, and Aq is taken to be zero, then we should 
get the classical result: in accordance with the "sum 
rule" the expected enhancement factor would be ~ 1 if 
LOc Al, and gc ~ Al/Ah if Wc ~ A^. But Aq is finite. 



(74) and we get 



5c [qm] 



1 



A:, V Ao 



X 5c [cl] 



(80) 



which is analogous to "weak localization corrections" to 
the mesoscopic conductance of closed rings [36j . 

Evaluation of gg.— The typical value of the elements, 
unlike the average value, is dominated by the majority 
of small elements. In order to calculates gs as defined in 
Eg. ([5^ . we have to bridge between the FOPT and the 
semiclassical analysis. To do it in a mathematically rig- 
orous way seems to be impossible. We therefore extend 
standard phenomenology and test it against numerical 
results. The basic idea is that FOPT cannot be trusted 
globally once levels are mixed non-perturbatively, but 
still it can be used in a restricted way. The analogy here 
is with Wigner's Lorentzian whose tails are given cor- 
rectly by FOPT, in-spite of the non-perturbative mixing 
of levels. See discussion of this issue in (33| . 

It is natural to expect FOPT to hold as an estimate for 
the majority of small elements as long as it does not ex- 
ceed the semiclassical estimate. If we take a band match- 
ing cutoff ujc ~ Apt; and calculate the ratio of the "area" 
under Eq.((75|) to the "area" under Eq.(|76|) we get: 



9s 



6-4a 



(81) 



Note that with a = 1 it follows that gs oc [u/K)^. In 
our numerics we fix uJc as the first minimum of Ca{^j->) 
implying ujc ^ Al, and consequently gc ^ 1, and g ^ g^. 
Our numerics fits well to 5 oc u'^ /h, indicating that the 
effective a is somewhat larger than unity. 

At this point one should appreciate how the contra- 
dicting FOPT and semiclassical results reconcile. The 
former apply to the majority of elements while the lat- 
ter apply to the algebraic average which is dominated by 



9 



relatively rare elements. The WQC regime where this 
picture is valid is further discussed in Sec. ()VII|) . 

For completeness one should be aware that the typical 
(median) value of the elements in X provides an under- 
estimate for the resister network average ((X))^. The 
reason is very simple: even if the matrix is very sparse 
(s 1) a network becomes percolating if the bandwidth 
is large enough. An RMT perspective [3l|, that uses a 
generalized variable-range-hopping approach, implies the 
following prescription; 



g I— max 1 1 , g exp y^— In 6 In g | 



(82) 



Here b = uJc/A.Q is the dimensionless bandwidth. This 
prescription allows to "correct" the result that has been 
deduced for g on the basis of a typical value estimate of 
the matrix elements. It is required if b is large. 



VII. THE WQC REGIME 

Quantum mechanics introduces in the billiard problem 
an additional frequency scale Aq that corresponds to the 
mean level spacing. We can associate with it the Heisen- 
berg time t^ = 27r/Ao. It is also possible to define the 
Ehernfest time which is required for the exponential 
instability to show up in the quantum dynamics. One 
can write 



= [log{i/h)]t^ 



(83) 
(84) 



where d—2. The traditional condition for "quantum 
chaos" is t^ ^ in , but if we neglect the log factor it is sim- 
ply in ^ ^H- This can be rewritten as ^ Aq, which 
we call the frequency domain version of the quantum 
chaos condition. Optionally one may write a parametric 
version of the quantum chaos condition, namely u ^ Uf,, 
where 



Ub 



[de Broglic deformation] 



(85) 



Note that it is the same as the QCC requirement of 
Eq. ((74)) . Namely, the frequency domain version of this 
condition implies that it should be possible to resolve the 
zero frequency peak of C(w) as in FiglHl while the para- 
metric version means that a de Broglie wavelength defor- 
mation of the boundary is required to achieve "Quantum 
chaos" . 

In practice we witness a WQC regime instead of hard 
chaos. We observe in the upper panel of Fig.[9] that gs 
is significantly smaller than unity, even for very small 
values of h for which u > Ub is definitely satisfied. For 
completeness we show in the lower plot additional data 
points in the regime u < Ub where this breakdown of QCC 
is not a big surprise. We conclude that QCC for u> Ub 
is restricted to C(aj), and does not imply Hard quan- 
tum chaos (HQC), but only WQC. In the WQC regime 



Cs{r) Ca{r) and consequently ^ 1, indicating spar- 
sity 

This emergence of the WQC regime can be explained 
by extrapolating FOPT considerations. If a wall of a bil- 
liard is deformed, the levels are mixed. FOPT is valid 
provided \Unm\ < Aq. This condition determines a para- 
metric scale Uc- If the unperturbed billiard were chaotic, 
the variation required for level mixing would be (37j 

Uc ~ h/ik^Lf^ = h^l'^ [not applicable] (86) 

This expression assumes that the eigenstates look like 
random waves. In the Wigner regime {uc < u < Ub) there 
is a Lorentzian mixing of the levels and accordingly 



# mixed levels k, [u/ucY 



[not applicable] (87) 



But our unperturbed (rectangular) billiard is not chaotic, 
the unperturbed levels of the non-deformed billiards are 
not like random waves. Therefore, the mixing of the lev- 
els is non-uniform. FigIS] illustrates the mixing vs u. 

By inspection of the Un^n .m^m matrix elements one 
observes that the dominant matrix elements that are re- 
sponsible for the mixing are those with large n^: but 
small I \. Accordingly, within the energy shell 



E„ 



E, the levels that are mixed first are those with 



maximal n^, while those with minimal are mixed last. 
The mixing threshold for the former is 



h/ik^L) = 



while for the latter one finds u^ ^ Tf' , which is much 
larger than Ub = h^. In our numerics g « u'^/h, implying 
that the WQC-HQC crossover is at 



(89) 



and not at Ub = h. Accordingly, the WQC regime ex- 
tends well beyond the traditional boundary of the Wigner 
regime, and in any case it is well beyond the FOPT bor- 
der Uc- 

WQC in broader perspective . In a broader 
perspective the term WQC is possibly appropriate also 
to system with zero Lyapunov exponent [t^=0Q), e.g. 
the triangular billiard [1^, and pseudointegrable bil- 
liards [iOl, and to systems with a classical mixed phase 
space. But in the present study we wanted to consider 
a globally chaotic system, under semiclassical circum- 
stances such that Ar is quantum mechanically resolved 
and QCC is naively expected. In this context there are of 
course other interesting aspects, such as bouncing related 
corrections to Weyl's law [2l| , and non- universal spectral 
statistics issues (see below), while our interest was with 
regard to the semi-linear response characteristics of the 
system. 

Spectral statistics in the WQC regime . The 

spectral statistics in the WQC regime has been studied 
in |22| concerning nearly circular stadium billiard, and 
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in |23j concerning circular billiards with a rough bound- 
ary. The model that we analyze is not identical, but can 
be regarded as a variation on the same theme. In FigfTUl 
we display some results for the level spacing statistics 
P{S), where the statistics is over Sn = {En+i — i?„)/Ao. 
It can be fitted to the cumulative Brody distribution 



FiS) = 1 



-bS" 



r 



q + 2 

g + 1 



9+1 



(90) 



which interpolates between the Poisson distribution 
{q=0) and with the Wigner surmise (q=l). This cumula- 
tive distributions can be transformed into linear functions 
T{x) = ln[-ln(l - F{c^))] with respect to the variable 
X = In(S'), and the fitting to our data gives g=0.38. 

Let us remind very briefly how the WQC border is 
determined in this context. It is convenient to describe 
the dynamics using a Poincare map, which relates the 
angle 6r of successive collisions (r = 1, 2, 3, • ■ •) with the 
piston. One observes that due to the accumulated effect 
of collisions with the deformed boundary, there is a slow 
diffusion of the angle with coefficient 



Dg - 

Accordingly, the classical crgodic time is 

T„ ^ 1/Dg ^ 



(91) 



(92) 



of a cold atoms that are trapped in an optical billiard. 
We shall identify it as the dimensionless absorption coef- 
ficient, and we shall inquire the feasibility of witnessing 
the quantum suppression factor which is related to the 
connectivity of the induced Fermi-Golden-Rule (PGR) 
transitions. 

LRT.— In linear response theory one has to know the 
following information in order to calculate the EAR: (i) 
The temperature T of the preparation; (ii) The spec- 
tral fluctuations (7(0;) of the system; (iii) The spectral 
content S{lu) of the driving ; Let us elaborate on the lat- 
ter. The RMS value of the vibrating wall velocity can be 
written as 



/[RMS] 



rA 



(95) 



where A is the amplitude of the wall movement. The 
power spectrum of f{t) has a spectral support ujc- To be 
specific let us assume that 



= / TJ— exp 



(96) 



The wall vibrations induce diffusion in energy space. 
Within LRT the diffusion coefficient is given by the Kubo 
formula, which in the following version can be regarded 
as an Einstein fluctuation-dissipation relation: 



and the quantum brcaktime due to a dynamical localiza- 
tion effect is 



Th - Dg/fv" ^ {u/hf 



(93) 



The border of the WQC regime is defined by the condi- 
tion T/i < leading to Eq. (|89|) . However, we would not 
like to over-emphasize this consistency because it is not 
a-priori clear that spectral-statistics and sparsity related 
characteristics always coincide. 

Intensity statistics in the WQC regime . WQC 

is also reflected in the intensity statistics of the wave- 
functions. If we had HQC we would expect Porter- 
Thomas (Gaussian) statistics and random wave corre- 
lations. The wavefunctions that we find do not look like 
random waves. In Fig.[TT] we show the statistics of the 
integrated intensity: 



|^(")(y)pdy = - 



1 



2En 



(94) 



Note that the total intensity, which is obtained by in- 
tegrating along the whole boundary with proper weight, 
gives unity, corresponding to the normalization of the 
wavefunction. 



VIII. THE HEATING RATE PROBLEM 

In this section we would like to discuss the physical sig- 
nificance of g with regard to the response characteristics 



D 



C{uj)S{uj)duj 



(97) 



The EAR per particle for strongly chaotic dynamics, as- 
suming that correlations between collisions can be ne- 
glected, is given by the wall formula 043 • Here we use 
the 2D version 



T 2T 



Sir Lx 



f = Go/2 



(98) 



Regarding the ballistic period as the time unit, and T as 
the energy unit, the dimensionless EAR is 



E 

tK 



37r2 \ A, 



2 .^x2 

L 



(99) 



In the quantum context the level spacing Ag sets the 
natural units for both energy and time measurements. 
Accordingly, we calculate the dimensionless quantity 



D 

A3 

^0 



37r2 



A, 
Ao 



(100) 



FGR.- The LRT formula Eq. (|97l) can be obtained 
from a classical derivation, say using a kinetics Boltz- 
mann picture, that does not assume applicability of the 
FGR picture. The same formula is obtained from FGR 
but with reservations that we illuminate in the next para- 
graph. It is therefore important to figure out the bor- 
der between the quantum FGR regime and the classical 
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Boltzmann regime. The strict FGR condition states that 
the near-neighbor transitions between levels should have 
a rate < Wc- Taking into account that the diffusion 
coefficient can be written as D ^ be x wq x Aq, where 
be = Wc/Ao, it follows that the strict FGR condition can 
be written as 



D 

^0 



< 



Ao 



power 



(101) 



with power= 2. But to witness FGR physics we can allow 
non-perturbative mixing on microscopic energy scales. 
The more careful analysis of [i^ leads to the same con- 
dition but with power= 3. 

SLRT.— It has been illuminated in a series of publi- 
cations (29l - [3l| that in the FGR regime one should refer 
in general to semilinear response theory. SLRT applies 
to circumstances in which the environmental relaxation 
is weak compared with the /(t)-induced transitions. In 
such circumstances the connectivity of the transitions 
from level to level is important, and the diffusion coeffi- 
cient is obtained via a resistor network calculation. Let 
us give a more precise quantitative description of this 
latter statement. The absorption coefficient G is defined 
via Eq.dT]). This is strictly analogous to Joule law: here 
the heating is due the vibration of the piston, while in 
the Joule-Drude problem it is due the oscillation of an 
electric field. The calculation of G can be done either 
within the framework of LRT using the Kubo formula 
(getting Glrt), or within the framework of SLRT (29l - [3l| 
using a resistor- network calculation (getting Gslrt). The 
correlations between collisions lead in the LRT case to a 
result that one can write as 



5c Go 



(102) 



where the expression for Go is implied by Eq. (|98p , and 
gc is defined as in Ea. (|33p . Similarly it is convenient to 
write the outcome of the SLRT analysis as follows: 



G, 



— gs Glrt — gs 9c Go — g Go 



(103) 



where g and g^ are defined as in Eq. p4p and Eq. 
If QCC considerations apply, then gs ^ 1 with small h 
dependent corrections as in Eq. (|80| . 

The resuhs of SLRT differ from those of LRT if the 
perturbation matrix is either sparse or textured, which 
is the case if we have WQC circumstances. The LRT and 
SLRT numerical results for gc and for g are displayed in 

Figia 



IX. EXPERIMENTAL MANIFESTATION OF 
QUANTUM ANOMALY 

With slight changes in notations which we find appro- 
priate for the experimental context, we summarize again 



the main parameters of the problem: 

oJl = ballistic frequency (104) 

clIr ~ Lyapunov ergodization rate (105) 

ujc = vibrations frequency span (106) 

ujQ ~ mean level spacing (107) 

The length scales are the linear dimension L, the 
de Broglie (thermal) wavelength Ae as determined by the 
temperature (calculated for E ^ T), and the radius of 
curvature of the walls R. The associated dimensionless 
parameters are: 

h = Xe/L = dimensionless Planck (108) 

u ~ L/R = deformation parameter (109) 

b ~ u/h ~ dimensionless bandwidth (110) 

a = A/ L ~ scaled vibration amplitude (111) 

Note that u determines the ratio cjr/wl, while h deter- 
mines the ratio wo/wl, hence b ~ uj-h/ujq. Our interest is 
in the non-trivial possibility /i < u <C 1, else uj^ cannot 
be resolved. 

The system.— Following [l^, [H, [l^ we consider 
^^Rb atoms (m = 1.4 x 10~^^kg), that are laser cooled to 
low temperature of T sa O.l/jK, such that the de Broglie 
wavelength is Ab = 1 /xw. The atoms are trapped in an 
optical billiard whose blue-detuned light walls confine the 
atoms by repulsive optical dipole potential. The motion 
of the atoms is limited to the billiard plane by a strong 
perpendicular optical standing wave. Assuming that the 
linear size of the billiard is L = 10 //to, the dimensionless 
Planck is h= 0.1 leading to Wl/wo = 30. Note that 

[27r]wE/(2L) == 220 Hz (112) 
wo= [2^]ft2^„,/(mL2) =7.5 Hz (113) 

where the [2tt] should be omitted for Hz units. Assum- 
ing 10% deformation the dimensionless bandwidth can 
be timed as & ^ 10. 

By modulating the laser intensity, one of the billiard 
walls can be noisily vibrated. We assume that the driv- 
ing is band- matched, i.e. coc ~ Wr. These are roughly 
the same parameters as in our analysis, for which we ex- 
pect gs ~ 0.1. 

The SLRT anomaly.— The common-wisdom expec- 
tation is that if QCC applies with regard to C{uj), then 
from Eq. (j97| we should get for the absorption coeffi- 
cient roughly the same result classically and quantum 
mechanically. SLRT challenges this expectation. It ap- 
plies to circumstances in which the environmental re- 
laxation is weak compared with the /(t)-induced tran- 
sitions. In such circumstances the connectivity of the 
transitions from level to level is important, and the LRT 
result should be multiplied by gs- 

In order to witness the SLRT anomaly, the driving 
amplitude A should be large enough so as to have a mea- 
surable heating effect, but small enough such that the 



12 



FGR condition is not violated. Disregarding prefactors 
of order unity it follows from Eq. ([99| and Eq. (|101|) that 
the requirements are 

> 10-3 (114) 
xa^ < (115) 

The first condition is bases on the assumption that it 
is possible to hold the atoms for a duration of ^ 1000 
bounces. Accordingly, there is a range where both condi- 
tions are satisfied, and there the SLRT anomaly should 
be observed, provided environmental relaxation effects 
can be neglected. 

It is worth noting that our theory for G is called SLRT 
because on the one hand S{ijj) h- c5(cj) leads to G i— cG, 
but on the other hand S{uj) H> 51(0;) + S2{oj) does not 
lead to G M> Gi + G2 . This semi-linearity can be tested 
in an experiment in order to distinguish it from linear 
response. 

X. BALLISTIC VERSUS DIFFUSIVE 
SCATTERING 

The EAR due to low frequency driving is determined 
by the couplings |i^„m|^ between nearby levels. Let us 
see how conflicting expectations with respect to its de- 
pendence on u reconcile by the analysis that we have 
introduced. For a small deformation FOPT implies that 
the couplings are oc it^, and hence 

E (xu'^ [FOPT expectation] (116) 

As u becomes larger the common expectation, based on 
Wigner theory, is to have Lorentzian mixing, leading to 
couplings oc 1/m^, and hence one expects 

Eccl/u^ [Wigner expectation] (117) 

In the formally equivalent problem of a conductance cal- 
culation this "Joule law" implies that the conductance is 
G oc 1/m^, where u represents the strength of the disor- 
dered potential. For the purpose of derivation, instead of 
using the FGR or Wigner picture, one can use the equiv- 
alent Drude picture, where the Born mean free path is 
i oc \/v? . On the other hand QCC considerations, based 
on Eg. ((27)) and using Ea. ((5T|) . imply that the couplings 
should be oc 1/u, and hence one expects 

Eoil/u [QCC expectation] (118) 

We therefore encounter here 3 conflicting expectations 
for the dependence of the EAR on the deformation pa- 
rameter. The analysis that we have presented resolves 
the conflict. Let us emphasize the main insights. 

Ballistic scattering.— We have assumed a smooth 
deformation: the worst case was a ~ 1, but more gener- 
ally we might have softer deformations with a > 1. Con- 
sequently, the mixing is not uniform: there are levels that 



are not mixed even if the perturbation is strong enough 
to mix some other levels. This leads to an interesting 
co-existence of Semiclassical theory and FOPT. Namely, 
we observe that the ((|F„,„|^))a agrees with Semiclassics, 
while {{\Fnm\'^)) s is given essentially by FOPT. The stan- 
dard Wigner theory does not apply, and the EAR would 
is oc or oc 1/u depending on whether LRT or SLRT 
applies: As the driving strength is increased we expect a 
crossover from LRT to SLRT. 

Diffusive scattering.— If the deformation profile 
Duiy) is erratic on sub Ab scale, then U is somewhat 
similar to the white disorder that has been analyzed in 
Rcf. [31I [ll] . Under such circumstances all the matrix el- 
ements of Unm are comparable. Consequently, one would 
observe Lorentzian mixing oc u^. Therefore C{uj) would 
have a Lorentzian peak of width oc u^, which differs from 
the semiclassical peak oc u. Furthermore, taking into 
account that the area under the central peak of C (r) re- 
mains the same irrespective of u, one deduces that 

{{\Fnm\'^))a/s « \ [Wigucr mixing] (119) 

and hence very different from both the FOPT predic- 
tion oc u^, and from the semiclassical expectation oc u. 
In other words - for diffusive scattering, unlike ballistic 
scattering, QCC does not apply. If U were like "white 
disorder" the quantum dynamics would be characterized 
by the Born mean free path, which is very different from 
the classical mean free path. 



XI. SUMMARY 

It is important to realize that we are studying in this 
work a driven chaotic system, and not a driven integrable 
system. Remarkable examples for driven integrable sys- 
tems arc the kicked rotator [sol - li^ and the vibrating 
elliptical billiard [i^l . In the absence of driving such sys- 
tems are integrable, while in the presence of driving a 
mixed phase space emerges. This is not what we call here 
weak chaos. Rather our focus is on completely chaotic 
systems that have a very small Lyapunov exponent com- 
pared with the ballistic scale. 

Weakly chaotic systems do not fit the common RMT 
framework. The Hamiltonian matrix of such a driven sys- 
tem does not look like one that is taken from a Gaussian 
ensemble, but rather it is very sparse. One can char- 
acterize this sparsity by parameters s and g that reflect 
the percentage of large elements, and their connectivity, 
respectively. For g we have used a resistor network calcu- 
lation that has direct relation to the semi-linear response 
characteristics of the system. 

We have highlighted that weakly chaotic systems pos- 
sess a distinct WQC regime, much wider than originally 
expected, where semiclassics and Wigner-type mixing co- 
exist. Then we discussed the implications of this obser- 
vation with regard to the theory of response. 
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The heating of particles in a box with vibrating walls 
is a prototype problem for exploring the limitations of 
linear response theory and the quantum-to-classical cor- 
respondence principle. In the experimental arena this 
topic arises in the theory of nuclear friction , and 
in the studies of cold atoms that arc trapped in op- 
tical billiards [l^-O- Mathematically it is related to 
the analysis of mesoscopic conductance of ballistic rings 
[38{ . In typical circumstances the classical analysis pre- 
dicts an absorption coefficient that is determined by the 
Kubo formula [Ulii-ill, leading to the "Wall formula" 
in the nuclear context, or to the analogous "Drude for- 
mula" in the mesoscopic context. The question arises 
[13, US [13, IsoI - IsgI I are there circumstance in which 
the quantum theory leads to a novel result that docs not 
resemble the semiclassical prediction. 

The low frequency driving that we assume is stochastic, 
rather than periodic. This looks to us realistic, reflecting 
the physics of cold atoms that are trapped in optical bil- 
liards with vibrating walls. It is also theoretically conve- 
nient, because we can use the Fermi-Golden-Rule picture. 
If one is interested in periodic driving of strictly isolated 
system, then there are additional important questions 
with regard to dynamical localization [39l - l42l . [57| , that 
can be handled e.g. within the framework of the Floquet 
theory approach. 

We predict that the EAR of a weakly chaotic system 
in the WQC regime would exhibit an SLRT anomaly: 
An LRT to SLRT crossover is expected as the intensity 
of the driving is increased; the linearity with respect to 
the intensity of the source is maintained but with a dif- 
ferent (smaller) coefficient; while the linearity with re- 
spect to the addition of independent sources is lost. This 
anomaly reflects that the absorption process in the meso- 
scopic regime might resemble a percolation process due 
to the sparsity of the perturbation matrix. In systems 
with diffusive scattering, that are in the focus of stan- 
dard condense matter textbooks, such an effect could not 
arise. 
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which corresponds to Eq. (j96|) . For characterization pur- 
pose we assume a band-matching weight function, mean- 
ing that be is chosen as the natural bandwidth of the 
matrix, corresponding to A^. The algebraic average is 
defined in the standard way: 



{{X))a 



-E 



Fiji — Tfi) X„ 



(A2) 



where TV is the size of the matrix, which is assumed to be 
very large. The algebraic average is a linear operation, 
meaning that 

((AX)) = A((X)) (A3) 
{{X + Y)) = {{X)) + {{Y)) (A4) 

There are different type of "averages" in the literature, 
such as the harmonic average, geometric average, and we 
can also include the median in the same list. All these 
"averages" are semi-linear operations because only the 
{{XX)) = X{{X)) property is satisfied for them. Irre- 
spective of the semi-linearity issue any type of average 
should satisfy the following requirement: if all the ele- 
ments equal to the same number, then also the average 
should equal the same number. 

In this paper we highlight a new type of average that 
we call a resistor-newtwork average. The defining pre- 
scription for its calculation is simple: given Xnm we as- 
sociate with it a resistor network Gnm via Eq.pi|). and 
define {{X))s as its inverse resistivity. 

There are a few cases where an analytical expression 
is available for the inverse resistivity G of a network 
Gnm- If only near neighbor nodes are connected, allow- 
ing Gn,n+i = g-n to bc different from each other, then 
"addition in series" implies that the inverse resistivity 
calculated for a chain of length N is 



G 



-E- 



(A5) 



If Gnm — gn-rn IS a function of the distance between the 
nodes n and m then it is a nice exercise to prove that 
"addition in parallel" implies 



G = 



OO 



(A6) 



Appendix A: The resistor-newtwork average 

We use the notation {{X)) in order to indicate the 
average value of its in-band elements. First we would like 
to define the standard algebraic average. It is essential 
to introduce a weight function that defines the band of 
interest. In the physical context this function reflects the 
spectral content of the driving sources. In practice we 
use rectangular or exponential weight function, say 



Note that in the latter case the resistor network average 
coincides with the algebraic average. In order to have a 
different result the diagonals of the matrix should be non- 
uniform, which is the case for sparse or textured matrices. 

In general an analytical formula for G is not available, 
and we have to apply a numerical procedure. For this 
purpose we imagine that each node n is connected to a 
current source /„. The Kirchhoff equations for the volt- 
ages are 



F{r) 



1 

2bc 



-\7i — 'rn\/bc 



(Al) 



(A7) 
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This set of equation can be written in a matrix form: 

GV = I (A8) 

where the so-called discrete Laplacian matrix of the net- 
work is defined as 



(A9) 



This matrix has an eigenvalue zero which is associated 
with a uniform voltage eigenvector. Therefore, it has 
a pseudo-inverse rather than an inverse, and the Kirch- 
hoff equation has a solution if and only if = 0. In 

order to find the resistance between nodes = and 
'T'our = N , we set /q = 1 and /^v = — 1 and /„ = other- 
wise, and solve for Vq and Vn- The inverse resistivity is 
G=[{Vo-VN)/Nr\ 



Appendix B: Intensity of fluctuations - optional 
derivations 



In this appendix we clarify the low frequency behavior 
of C{u}) using two optional approaches. We assume that 
7 7o is roughly a constant, so there is a well defined 
correlation time 



1 

iR = - 

7 

Time domain approach. Observe that 



Fe{t)dt 




{qe/Te)t, 
lre)t, 



t < 

t > tR 



(Bl) 



(B2) 



From Eq.dlS]) it follows that (F(O)F(t)) looks as follows: 
at i = it contains a self-correlation delta peak; within 
t -C ia it is the theta averaged comb of delta peaks due to 
bouncing; For t 3> it flattens and reflects the squared 
average value of F . Accordingly, the short time average 
and the long time average values of {F{0)F{t)) are cq 
and Coo of Eg. ([38) ) and Eq. (p9|) . Consequently, the "area" 
under the correlation function is 



C{uj=0) w (co-Coo) X 2iR 
in agreement with Eq.(|51|). 



(B3) 



Number variance approach. It is instructive to 
deduce C(lu ~ 0) using over-simplifled derivation via the 
number variance approach, as in the analysis of spec- 
tral rigidity [s^. This over-simplifled approach treats 
the spikes as having equal size (below q = 1). The vari- 
ance in the number of collisions during the time t is given 
by the expression 

Var(iVm) = — / —^sm^(nojt)duj (B4) 
Jo ^ 

Consequently, 

C{uj = 0) = ^^^(^i^)) ^ diffusfon in counting (B5) 



Assuming that the step in this random walk process is of 
duration ^r , the diffusion coefficient is 



= 0) = -Var ( - 



ivar f i 
7 \T 



(B6) 



Which leads upon restoration of q to Eq. ([5T 



Appendix C: The Fnm matrix for zero deformation 

Here we calculate the large scale sparsity po, and the 
average value of |F„mp, in the case of a rectangular box. 
It is tempting to identify pq as " s" , but in fact the latter 
is ill defined because it refers to the sparsity of the in- 
band elements, while for u = the bandwidth is zero. 

We consider the matrix elements that reside in- 
side an energy window of width 6E. The levels 
(jixjUy) within this window belong to the energy shell 
E < En^,ny <E + 6E. We define the "radius" of this 
shell as = (2mi?)^/^. For a given riy section the width 
of the shell is denoted as Sn^., and in wavenumber units 
it is given by the expression 



1,2 



kpSkiT 



1.2 



h.2 



(Cl) 



The total number of levels within this window can be 
calculated in a complicated way as 



Sn^dny 



LxLy 



m5E _ SE 

~ IT' 

^0 



1,2 _ h.2 

l^y 



(C2) 
(C3) 



Similarly we can calculate the number of coupled levels, 
and hence the large scale sparsity: 



Po 



1 



V Jo 



1,2 _ U2 *!/ 



(C4) 



Trk^Ly 



In 



dk 



The FOPT perturbed matrix is sparse and textured. 
Its non-zero elements are of size k^/mLx as implied by 
Eq. (|6T|) . The algebraic average of the elements is given 
by 



{{\Fnm I ))oo — 



Sri^dny 



4 
ttL. 



y Jo 



T2 u'^^y 



' k 


X 






2 


kl 





(C5) 
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kl 



which is the same result Eq. 
estimate. 



as in the semiclassical 
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FIG. 1: (Color online) Sketch of the billiard system of Eg . (fT3)l . 
The unperturbed billiard is a rectangle of size x Ly. The 
deformation U, due to the radius of curvature R of the left 
wall, is characterized by the parameter u = Ly/ R. In order 
to break the mirror symmetry the center of the curved wall 
is shifted upwards a vertical distance e. The time dependent 
perturbation is due to the displacement f(t) of the right wall. 




FIG. 2: (Color online) The driving induces transitions be- 
tween levels En of a closed system, leading to diffusion in 
energy space, and hence an associated heating. The diffusion 
coefScient De can be calculated using a resistor network anal- 
ogy. Connected sequences of transitions are essential in order 
to have a non-vanishing result, as in the theory of percolation. 
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FIG. 3: (Color online) An image of the eigenstate i5„~13618 
for the billiard of Fig[l]with Lx = 1.5, and Ly=1.0, and R—8, 
and e = 0.1. In the numerics the units are chosen such that 
ftpianck^l and the mass is m=l/2. 
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FIG. 4: (Color online) Using a truncated matrix represen- 
tation of the deformed billiard Ho -I- 17 in the unperturbed 
basis n = {nx,ny) of the non-deformed rectangular billiard 
Ho, we find a representative eigenstate no of the former and 
plot \{Efi\E„g)\'^ versus the running index n of the ordered 
energies. The participation number (FN) of an eigenstate in 
this basis reflects how many energy levels were mixed due 
to the deformation u. Having PN>1 indicates that our data 
is beyond the FOPT regime. For the displayed eigenstate 
PN ~ 5.77, while its energy width is 19 levels. The billiard 
parameters are R = 8 and 1/ft « 27.15. 




FIG. 5: (Color online) The participation number (PN) of the 
eigenstates as a function of u (panel a) and h (panel b). The 
parameter h characterizes an energy window that contains 
100 eigenstates. The method of calculation is as explained 
in Fig.|4l The average energy width in units of mean level 
spacing is represented by red solid line. Having average width 
larger than the average PN is an indication for sparsity. 
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FIG. 6: (Color online) Image of the perturbation matrix X — 
{\F„mf} for the billiard of FigEl within the energy window 
3500 < E„ < 4000. This matrix is sparse. More generally 
it might have some texture. The latter term applies if the 
arrangement of the large elements is characterized by some 
pattern. 
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FIG. 8: (Color online) The band profile of the perturbation 
matrix for the billiard of Fig.[3] where R = 8. (a) The al- 
gebraic average and median along the diagonals of the X„,n 
matrix versus uj = {En— Em). The vertical axis is normalized 
with respect to Coo, while the horizontal axis is uj/v-e. The 
classical power spectrum is presented to demonstrate the ap- 
plicability of the semiclassical relation Eq. (|27|I . The red line 
is the analytical expression that applies to zero deformation. 
The dotted vertical line is the frequency 1/ti^ and the dashed 
one is l/tn. (b) Zoom of the ui <C 1/ti^ region. For sake 
of comparison we display results also for R = 2. The vertical 
lines indicate the mean level spacing. The dashed red curves 
are based on EQ. (f52|) . 



FIG. 7: (Color online) Histogram of the values of the ele- 
ments Xnm that reside in the central band of the matrix. 
The analysis is done for the billiard of Fig.[3] where R — 8. 
The statistics includes all the elements in the energy windows 
100 < -B < 4000 (EWl), and 10000 <E < 14000 (EW2). For 
sake of comparison we display results also for R = 1. 
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FIG. 9: (Color online) SLRT vs LRT. The scaled absorp- 
tion coefficient Qc (LRT) and g = QsQc (SLRT) versus the 
dimensionless 1/h (upper panel), and versus the dimension- 
less deformation parameter u = L/ R (lower panel) . Note that 
t; = 1 is the prediction of the "Wall formula" , while the line is 
based on the classical analysis. In the upper panel the anal- 
ysis has been done for the billiard of Fig{Tl The calculation 
of each point has been carried out on a 100 x 100 sub-matrix 
of X centered around the h implied energy E. The "un- 
textured" data points are calculated for an artificial random 
matrices with the same bandprofile and sparsity (but no tex- 
ture). The complementary lower panel is oriented to show the 
small u dependence. The analysis is based on a truncated ma- 
trix representation of Ho + U , within an energy window that 
corresponds to 1/ft ~ 9. Due to the truncation there is some 
quantitative inaccuracy with regard to the larger g values. 



FIG. 10; (Color online) (a) Cumulative histogram for the level 
spacing distribution P{S) with fitting to Brody distribution 
(g=0.38), and contrasted with Poisson distribution (g=0) and 
the Wigner surmise {q=l). (b) The Brody parameter is de- 
termined via the slope of T{x) as explained in the text. 
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FIG. 11: (Color online) Wavefunction intensity statistics: 
The cumulative distribution of the integrated intensity /„ of 
Eg . (|94p is presented. As R becomes larger the distribution fur- 
ther deviates from Gaussian statistics. The statistics includes 
all the eigenfunctions in the energy windows 100 < E < 4000 
(EWl), and 10000 <E < 14000 (EW2). 



